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Abstract 


A production 3-D elliptic flow code has been used to calcu- 
late non-reacting and reacting flow fields in an experimental 
mixing section relevant to a rich-burn/quick-mix/Iean-burn 
(RQL) combustion system. A number of test cases have been 
run to assess the effects of the variation in the number of 
orifices, mass flow ratio and rich-zone equivalence ratio on 
the flow field and mixing rates. The calculated normalized 
temperature profiles for the non-reacting flow field agree 
qualitatively well with the normalized conserved variable 
isopleths for the reacting flow field indicating that 
non-reacting mixing experiments are appropriate for screen- 
ing and ranking potential rapid mixing concepts. For a given 
set of jet momentum-flux ratio, mass flow ratio, and density 
ratio (J, MR, and DR), the reacting flow calculations show a 
reduced level of mixing compared to the non-reacting cases. 
In addition, the rich-zone equivalence ratio has noticeable 
effect on the mixing flow characteristics for reacting flows. 


Nomenclature 


A area 

A m duct crossectional area, also A (ot , m 2 

ACd also Aj, effective orifice area m 2 
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area ratio (jet/mainstream) = Aj/A m = ACd/A m 
area weighted overall T or <j) deviation from 
equilibrium, Eq 3 

area weighted T or <(> deviation above equilib- 
rium, Eq4 

area weighted T or <j> deviation below equilibrium, 
Eq 5 

area determined half width of the distribution 
function 

constant of proportionality between V J and S/R. 


Eq-9 


diameter of the orifice 

density ratio (jet/mainstream) 

total pressure loss across the mixing wall, % 

non-dimensional temperature and/or equivalence 

ratio, Eq 2 

fuel to air ratio 

radial vector direction 

momentum-flux ratio (jet/mainstream) = M 2 /DR, 


also (MR) 2 /[(DR)(A/A m ) 2 ] 

tangential vector direction 

overall mixer mass flowrate, kg/sec = mj et 

m main 

mass flux ratio (jet/mainstream) = DR Vje, 


+ 


mass flow weighted overall T or <|> deviation from 
equilibrium, Eq 6 

mass flow weighted T or $ deviation above 
equilibrium, Eq-7 

mass flow weighted T or <J> deviation below equi- 
librium, Eq-8 

mass flowrate ratio (jet/mainstream) 
optimum number of orifices / row, Eq-9 
total pressure, atm. 
radius of the mixing section, m 
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r radial distance from the centerline of the mixer, m 

S orifice spacing in the circumferential direction 

T temperature. K 

u local axial velocity, m/sec 

U ma j n approach mainstream axial velocity, m/sec 

VR velocity ratio (jet/mainstream) =V jct /U main 

V je , radial velocity ol the jet. m/sec 
w jet width, m. see Fig-3 

x axial distance from the leading edge of the orifice 

a parameter used in Eq-2 through 8. can be either 

temperature or equivalence ratio 
8 value of the cumulative volume fraction at the 

fequil . Fig-4 

p fluid density, kg/m 3 

<j) equivalence ratio ( f/a) tocal / ( f/a) stoj 


Subscripts 


equil equilibrium 

j jet 

m mainstream, also (main) 

RZ rich-zone 

stoi stoichiometric 


Introduction 


Advanced fuel-efficient commercial turbo propulsion gas 
turbines pose a number of design challenges for combustion 
systems including durability, aerothermal performance, wide 
operability range, and exhaust emissions. Conventional gas 
turbine combustion system employs a single stage combus- 
tion zone which has been optimized to produce low exhaust 
emissions of carbon monoxide, and unburned hydrocarbons 
at idle and near-idle operating conditions. Exhaust smoke 
below the threshold of plume visibility throughout the engine 
operation range has also been maintained without adversely 
affecting engine starting, stability, and relight requirement. 
Moderate reductions in NO x emissions (up to 30%) have 
been achieved in single-combustion-zone combustors by 
reducing unmixedness and combustion zone residence time. 
Further reductions in high-power NO x emissions (e.g. up to 
50% from state-of-the-art levels) would require design and 
development of two-stage combustion concepts including 
radially or axially staged combustion zones. These types of 
combustion systems are currently under development. 

In order to achieve more demanding (70 to 90%) NO x reduc- 
tion goals of the advanced turbine engines which will also 
have considerably higher overall pressure ratios and turbine 
rotor inlet temperatures, significant advances are needed in 
the design and development of combustion systems employ- 
ing multi-staging concepts including lean premix (LPP) and 
rich burn/quick mix/lean burn (RQL) (Shaw. 1991). The 
RQL concept studied in-depth by Rizk and Mongia (e.g. Rizk 


and Mongia. 1990; 1991; 1993) appears to have a number of 
advantages over other competing ultra-low NO x concepts 
provided one can significantly enhance the mixing during the 
transition from rich to lean-side combustion. Additional 
experimental and analytical research effort (e.g. Howe., et 
al.. 1991; Smith. Talpallikar. and Holdeman. 1991; Talpal- 
likaretal.. 1991: Vranos et al.. 1991:Bain. Smith. and Hold- 
eman. 1992. 1993: Liscinsky et al.. 1992: Zhu and Lai. 1992; 
Doerr and Hennecke. 1993: Liscinsky. True, and Holdeman. 
1993) is underway to study and identify the critical design ^ 
and flow parameters affecting the mixing effectiveness. 

An integrated analytical and experimental investigation . 
between the authors and the UCI Combustion Laboratory has 
been underway to quantify non-reacting mixing processes in 
a 80 mm diameter cylindrical test section, as summarized by 
Hatch et al.. 1992. Kroll et al.. 1993. and Oechsie. Mongia. 
and Holdeman. 1992. 1993. These studies have shown that 
the mixing section effectiveness is affected by the 
jet-to-mainstream momentum-flux ratio (J), mass flow ratio 
(MR), and orifice/slot design parameters including R/d, and 
S/d. aspect ratio and slot orientation. Here R. d. and S are the 
radius of the mixing section, orifice diameter, and circumfer- 
ential spacing between the orifice centers, respectively. The 
qualitative agreement between predictions and measurements 
is reasonably good, and the model calculations can be used to 
provide further insight especially for the flow variables that 
have not been measured. 

A logical question that could be raised is: are the 
non-reacting mixing investigations appropriate for defining 
an optimum mixer geometry. This paper summarizes the 
results of a parallel study to compare the mixing performance 
of an RQL mixing section in both reacting and non-reacting 
environments using constant main flow parameters such as 
jet-to-mainstream momentum-flux ratio (J). mass flow ratio 
(MR), and density ratio (DR) and establish the difference in 
the flow field structure. In this study, the mixing effective- 
ness of a jet in a crossflow is numerically investigated using a 
3-D analytical code. Eight different circular orifice configu- 
rations were analyzed in a cylindrical mixing section at J=25 
with both 8 and 12 orifices per row. The reacting flow 
configurations were calculated with mainstream incoming 
species mole fractions simulating reacting rich-zone equiva- 
lent ratios of 1.35 and 1.80. The procedure used to model the 
geometric configurations was similar to that used by Oechsie, 
Mongia. and Holdeman. 1992. 1993. 


Mathematical Model 


A production 3-D combustor code. COM-3D (Bruce. Mon- 
gia. and Reynolds. 1979) is used that solves the turbulent 
reacting flow transport equations using the SIMPLE algo- 
rithm of Patankar and Spalding (Patankar. 1980). This 
program simulates turbulence by the two-equation k-e model 
(Launder and Spalding. 1974). and combustion following 
vaporization is determined by a four-step chemical reaction 




model based on Arrhenius and modified eddy breakup con- 
cepts. The transport equations for all dependent variables are 
of the following form as shown in Eq- 1 : 

div[p r u^- (Mcff/Pr) grad(^)] = (1) 

where p r is the mixture density, u is the velocity. p eff is the 
effective turbulent viscosity. P r is the effective 
Prandtl/Schmidt number, and Sc is the source term for the 
variable q. The following variables are computed by 
COM-3D: 1) axial, radial, and swirl velocity components: 2) 
specific enthalpy and temperature; 3) turbulence kinetic 
energy and dissipation rate: 4) unburned fuel. CO. Hi. inter- 
mediate fuel, and composite fuel mass fractions: and 5) fuel 
spray trajectory and evaporation rate. 

The computational effort is significantly reduced by model- 
ing a sector of the mixing section comprising a single orifice. 
Therefore, t lie shape of the sector was dependent on the 
number of orifices equally spaced in the circumferential 
direction. It should be noted that 8 orifices/row yield a 
computational domain of a 45 degree sector, and similarly, 
12 orifices/row yield a 30 degree sector. Periodic boundary 
conditions were applied in the circumferential direction. 
No-slip and adiabatic boundary conditions were applied at 
the outer wall defining the inside wall of the mixing section. 
Zero-gradient boundary conditions were applied at the center 
axis. Axial gradients at the exit boundary condition were 
assumed zero. 


Geometric Configuration 

In this study, the mixing section was modeled as a constant 
diameter cylindrical duct with a single row of equally spaced 
orifices. The outer wall diameter is 80 mm and the axial 
length of the mixing section extended from x/R=-1.4 to 
x/R=6 where x is referenced from the leading edge of the 
orifice. Sufficient axial distance was provided both upstream 
and downstream of the orifice to minimize the impact of the 
inlet and exit boundary conditions on the calculated flow 
structure in the primary domain of interest which is 
0<x/R<l. The computational grid domain was typically 
discretized into 50.000 to 80,000 finite control volumes. The 
grid was typically denser near the orifice and near the outer 
wall to resolve the high velocity and temperature gradients 
resulting from the inlet of the crossflow jet. An orthogonal 
view of a typical grid arrangement is shown in Figure- 1. The 
grid is normally configured to allow smooth progressive 
volume change between adjacent control volumes to help 
speed up the convergence of the solution. 

The geometric configurations of the jet orifices is also shown 
in Figure- 1. A total of 8 circular hole configurations were 
analyzed. Also, it should be noted that the normalized tem- 
perature. normalized equivalence ratio, and velocity profile 
plots shown in this paper depict plane sections in a) the 


axial-radial plane through the center of the jet from x/R=-0.5 
to x/R=2.5, and b) the radial-tangential plane at x/R=l (one 
full mixer radius downstream from the leading edge of the 
jet). The axial, radial, and tangential directions are shown in 
Figure- 1. It should be noted that the flow direction is from 
left-to-right in the axial-radial planes and into the paper for 
the radial-tangential planes. 

The overall description of the eight configurations is given in 
Table- 1 and Figure-2. Cases 1. 3, 5. and 7 are for 
non-reacting flows and the corresponding reacting flow cases 
are denoted by 2. 4. 6. and 8. respectively. It is assumed that 
practical RQL mixing sections will have liner pressure drop 
and J approximately 3% of P mam and 25. respectively. The 
rich-zone equivalence ratio is assumed to lie between 1.35 
and 1.80 and lean-zone equivalence ratio is typically 0.45. It 
is also assumed that the non-reacting experiment is per- 
formed at atmospheric pressure whereas the reacting experi- 
ment might be run at high inlet temperature and pressure 
conditions. It is further assumed that both reacting and 
non-reacting experiments will be conducted with comparable 
values of J. mass flow ratio (MR), density ratio (DR), area 
ratio (Ar). and quick mix section liner pressure drop (DP/P). 

The eight cases that have been run are for fixed values of 
J=25: MR=1.83 and 2.67: corresponding to rich-zone equiva- 
lence ratio. 4)1*2=1.35 and 1.80. and other variables listed in 
Table-1. The blockage is defined as the circumferential 
projection of the orifice divided by the spacing between the 
orifice centers. Typically 150 to 200 finite control volume 
nodes were used to simulate the orifice. 

The reacting flow conditions, the operating pressure and 
temperature simulated realistic gas turbine operating condi- 
tions (P= 13.4 atm.. T= 978K). The mainstream flow and jet 
flow conditions are show in Table- 1. The incoming main- 
stream species mole fractions for (CO. COj. HiO. and Hi) 
were determined based on the given rich-zone equivalence 
ratio (see Table-1) at equilibrium at the prescribed operating 
conditions using JP-5 for the typical fuel properties. 
Although this assumption may not be totally realistic in an 
operating RQL rich section, it does provide a baseline for 
comparison. It is recognized that the rich-zone will produce 
a great amount of CO and unburned hydrocarbons which are 
directly admitted to the mixing region. This will probably 
increase the reaction in the mixing zone prompted by the 
mixture with the additional air. However, since the rich-zone 
performance greatly depends on the residence time and 
geometric parameters, and since its design is not within the 
scope of this study, the assumed rich-zone performance 
(achieving chemical equilibrium) is a reasonable initial 
condition towards analyzing the reaction in the mixing zone. 
The four specified species used in the inlet mainstream 
condition correspond to the four-step chemical reaction used 
in COM-3D. The species mole fractions, temperature, and 
velocity profiles were assumed to be uniform across the inlet 
cross-section of the mixing section. The jet-to-mainsteam 
mass flow ratio was determined by the prescribed lean zone 
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equivalence ratio of 0.45 (which is constant for all reacting 
flow configurations reported herein). The air jet flow was 
characterized by a radial, uniform flow across the orifice 
effective area. The assumption of uniform mass injec- 
tion/area is applied in the mathematical model in all the 
analyzed configurations. 

The non-reacting configuration models were setup for typical 
atmospheric temperature and pressure conditions (Tj et = 
300K. P=1 atm.) for which the mainstream inlet flowrate was 
varied to obtain the same J. MR. DR. orifice ACd. and DP/P 
as the comparable reacting How case. Note that in the actual 
optimization process, slight variations in the density ratio 
were obtained as shown in Table- 1. 

The turbulence kinetic energy of the mainstream and jet 
flows were 0.3% of the square of the mean velocities. The 
turbulence length scales of the mainstream How were 2% of 
the mixer diameter, and the turbulence length scale of the jet 
was of the order of the orifice diameter. The inlet conditions 
for all orifices in the mixing section were equal to create a 
symmetrical input condition about t lie circumferential direc- 
tion which was necessary for the sector analysis. 

A typical numerical solution took about 250 iterations for full 
convergence with overall mass flow residuals of 0.05% of the 
total mixing section mass flowrate. All solutions were 
obtained using the Cray Y-MP and a typically converged 
solution took about 1 to 1.5 hrs of CPU time. The reacting 
flow solutions consumed about twice the CPU time due to the 
additional species variables which are solved in the model. 


Results and Discussion 


The mixing performance for all configurations analyzed in 
this study were ultimately evaluated at x/R =1. It is however 
recognized that the mixing performance throughout the 
mixing section volume of interest (0<x/R<l) should also be 
considered since complex structures are present in the flow 
field especially near the entry of the jet. 

Two different methods were used to analyze the numerical 
results: 

l) the results were analyzed qualitatively by visual observa- 
tion of the temperature, equivalence ratio, and velocity field 
solutions. The temperature and equivalence ratio results are 
presented as normalized values with respect to the overall 
differential between the mainstream flow and the jet flow 
inlet temperatures and equivalence ratios consecutively. This 
normalized parameter (f) is defined in Eq-2 where a = tem- 
perature or equivalence ratio. 


The value of (0 varies from 0 to 1. where 0 is the value of the 
unmixed jet and 1 is the value of the mainstream flow. Note 
that f = 1-0. where 0 is as defined previously (Holdeman. 
1993) and used elsewhere also. It is also important to note 
that the definition of (f) also applied for the normalized 
equivalence ratio in the reacting flow configurations where 
f=l equals the rich-zone equivalence ratio and f=0 is the jet 
equivalence ratio. The definition of (f) however does not 
apply to the temperature distribution in the reacting flow 
solutions since the temperature in some areas of the mixing 
region rise above that of the incoming mainstream flow due 
to the chemical reaction of the rich-zone exit composition. 
The value of the equilibrium temperature (as applied to the 
non-reacting conditions) and equilibrium equivalence ratio 
(applied to the reacting conditions) based on the normalized 
definition parameter (0 in Eq-1 are also equivalent since both 
are conserved variables. Both equilibrium values were 
calculated in a purely adiabatic system at any location down- 
stream of the jet injection. 

In addition to the temperature and equivalence ratio distribu- 
tion plots, the velocity plots in the axial-radial plane through 
the orifice center are shown for all configurations for x/R 
from -0.5 to 2.5. In each plot, the streamlines corresponding 
to the developing jet in the mixing section are shown depict- 
ing the jet penetration and jet width. The values of the jet 
width (w) and mean jet penetration evaluated at x/R= 0.625. 
1 , and 2 are tabulated in Table-2 and the pictorial definition 
of both of these parameters is shown in Figure-3. The jet 
penetration at x/R=2 is shown only for the purpose to estab- 
lish if complete penetration is obtained at x/R=l. 

2) Three different statistical methods are reported in this 
paper: 

a) The performance of the mixing section at the 
radial-tangential planes at x/R=l was quantified by using 
area weighted planar deviation parameters. The smallest 
deviation indicates the best mixing configuration. Three 
different parameters (mixing deviations) AMIX. AHOT. and 
ACOLD are described in Eq-3 through Eq-5. Note that 
AMIX also equals the square root of the sum of the squared 
values of both AHOT and ACOLD. It is also important to 
note that the definitions in Eq-3 through Eq-5 do not correct 
the mixing non-uniformity for the bias introduced in the 
region of the section where the air is being injected through 
the orifice. This is accepted since the comparison between 
the reacting and non-reacting results is only performed at 
x/R=l which is downstream of the trailing edge of all ori- 
fices. The area-weighted non-uniformitv results are shown in 
Table-3. 
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b) Mass How weighted planar deviation parameters (also 
evaluated at x/R=l) as defined in Eq-6 through Eq-8 were 
also used to evaluate the mixing region. The MMIX value 
also equals the square root of the sum of the squared values 
of both MHOT and MCOLD. Note that these parameters are 
essentially the same as the area weighted parameters with the 
added density and velocity weighting terms. 


shown in this paper is the integrated histogram for (0<f<l). 
Therefore this value is generally independent of the bin size 
used in this analysis. Although not compared in this paper, 
the cumulative volume fraction for the interval f equil ± Af 
may be used to quantify the mixing uniformity within the 
tolerance Af. In this paper however, the cumulative volume 
fraction is only used to indicate whether the histograms 
indeed integrate to a value of 1 . as expected. The shape of the 
volume fraction histogram was also characterized by the 
definition of B(+) and B(-) similar to that used by (Oechsle. 
Mongia, and Holdeman. 1992). The value of B(+) is the 
"area determined" distribution half width above f equ ;i such 
that the integrated area under the histogram above feq U ,i is 1/2 
that of the overall histogram area above f equil . Tne same 
applies to B(-) for the area below f eqJlll . The definitions of 
B(+) and B(-) are shown pictorially in Figure-4. The best 
mixing results yield B(+) = B(-) = 0. meaning that the histo- 
gram distribution width effectively collapses towards the 
equilibrium value of (f) due to a homogeneous mixture of 
either temperature or equivalence ratio. 
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c) The flow field was also evaluated by performing a 
numerical volume integration throughout the mixing section 
of interest. The volume integration parameters allow the 
quantification of the entire (low field mixture which is more 
descriptive of the overall flow phenomena in the mixing 
section as compared to a planar deviation analysis as per- 
formed with the mass flow weighted and area weighted 
parameters. This volume integration was performed on the 
temperature or equivalence ratio distributions depending if 
the configuration analyzed was either non-reacting or react- 
ing respectively. The entire range of the normalized param- 
eter f (from 0 to 1 ) was sub-divided into 200 equal size bins 
and the volume of the computational control volumes corre- 
sponding to the value of (f) at a certain bin size (fj to fi+Af) 
was integrated. The integrated volume in each bin was 
normalized based on the entire analyzed mixer volume, thus 
obtaining the normalized volume fraction. The volume 
fraction for each bin was plotted in the ordinate with the 
corresponding value of f in the abscissa forming a histogram 
plot. In addition, the incremental partial volume per bin was 
integrated from f =0 to 1 and the cumulative volume fraction 
was obtained for all the 8 analyzed configurations. The 
cumulative volume fraction is also plotted versus (f) for all 
configurations. Note that the cumulative volume fraction 


The results in this paper are presented in the following man- 
ner: 

a) effect of reaction on mixing 

b) effect of number of orifices/row on mixing, and 

c) effect rich-zone equivalence ratio on mixing 


Effect of reaction on mixing 

The normalized temperature distribution results for all con- 
figurations for the axial-radial plane sections through the 
center of the jet are shown in Figure-5. The corresponding 
radial-tangential planes at x/R=l are shown in Figure-6. 
Significant differences are apparent in the temperature 
distribution for the same configuration with similar operating 
conditions for both reacting and non-reacting cases. Note 
that the approach flow normalized temperature in the react- 
ing flow cases is not the maximum value in the mixing 
section since the latter depends on the rich-zone equivalence 
ratio and is therefore not a conserved scalar. A lower 
approach mainstream normalized temperature value is 
expected with higher rich-zone equivalence ratio. These 
major differences are attributed to the fact that the reacting 
flow configurations are not only mixing hot and cold jets but 
also chemically reacting in areas permissible by the chemical 
kinetics of the species exiting the rich-zone of the combustor. 
It is worth noting that significant reaction occurs behind the 
orifice where the jet wake produces a low velocity recircula- 
tion between hot rich mixture exiting the rich-zone and 
relatively cold jets. This observation is shown for configura- 
tion # 2 and 4 (Figure-5b. 5d. 6b. and 6d) especially since 
both configurations have only 8 orifices. The configurations 
# 6 and 8 (Figure-5f. 5h. 6f. and 6h) however have 12 
holes/row thus producing a much smaller jet trailing wake 
and therefore significantly inhibiting this reaction in the 


mixing section. Note that the jet penetration tor 8 orifices is 
greater than for 12 orifices and this is discussed later in more 
detail. The non-reacting flow fields are similar to the results 
reported by Occhsle. Mongia. and Holdeman 1992. and 1993. 
Note that in Figurcs-5 and 6. f=l applies to the approach 
mainstream normalized temperature in the non-reacting flow 
configurations. Similarly, the stoichiometric temperature in 
the reacting flow configurations has a value of 1. 

The normalized temperature distribution (f) in the 
non-reacting flows compare well with the normalized equiva- 
lence ratio distribution for all the configurations in Figures-7 
and 8. These results seem to indicate that non-reacting 
temperature profiles can emulate the reacting flow mass flow 
distribution reasonably well. It is worth noting that the 
non-reacting jets appear to interact more at the center of the 
mixer at about x/R=0.5 as compared to the reacting flow jets, 
see (Figure-7a. 7b). This interaction is usually spotted by the 
upstream swirling flow produced at the location where 
opposing jets merge thereby inducing extra mixing between 
the mainstream and jet. This was not observed for the 12 
orifice/row configurations due to shallow jet penetration 
produced by the 12 orifice/row configurations. Note that in 
Figure-7 and 8. f=0 represents the jet normalized temperature 
(non-reacting flow), and f=l represents the mainstream 
normalized temperature (non-reacting flow). 

The planar statistical variances indicate very similar conclu- 
sions to the observed temperature and equivalence ratio 
distributions. The area weighted and mass flow weighted 
parameters in Table-3 indicate very similar mixing 
non-uniformities between reacting and non-reacting cases for 
the same J. MR and DR at x/R=l. These results are also 
shown graphically in Figure-9a and 9b. The results indicate 
that the non-reacting flow configurations mix better com- 
pared to the reacting flow counterparts. Both area weighted 
and mass flow weighted results indicate similar results as 
shown previously in Oechsle. Mongia. and Holdeman. 1993. 
The volume fraction histogram results for all the configura- 
tions are shown in Figures-9c. 10. and 1 1. The B(+) and B(-) 
half width parameters are also tabulated in Table-3 and both 
indicate somewhat similar conclusions to the area weighted 
and mass flow weighted parameters; however, the differences 
between the volume fraction histogram shapes for reacting 
vs. non-reacting flow configurations with similar geometry 
and operating conditions are plotted in Figure- 12. The 
results shown in Figure- 12 indicate that for all comparable 
configurations the cumulative difference of the volume 
fraction/bin amounts from 50% to about 70% of the overall 
mixer volume of interest see Figure- 12a. 12d. This observa- 
tion is of significant value since particular mixing flow 
configurations may compare well at x/R=l using both visu- 
ally and planar average methodologies at x/R=l but may be 
significantly different throughout the mixing section as 
shown by the volume integral results. This observation may 
have a significant effect on the future use of non-reacting 
flow configuration analysis when selecting a good mixing 
configuration for low NO x application. 


The penetration of the jet for both reacting and non-reacting 
configurations are tabulated in Table-2 and the results are 
shown in Figure- 13. The results indicate very little differ- 
ence in both the penetration of the jet core and the jet width 
at x/R of 0.625. 1. and 2. It is important to note that although 
the change in the jet width is small according to the results 
shown in Table-2, the temperature and equivalence ratio 
cores of the jet as shown in Figure-6 and 8 for the reacting 
cases are significantly different compared to the non-reacting 
cores (noting that J. MR. and DR in the comparison are 
constant). Both temperature and equivalence ratio cores of 
the jet are significantly preserved throughout the flowfield up 
to x/R=l for the reacting flow calculations only. 


Effect of the number of orifices/row on mixing 

The effect of the number of orificcs/row on mixing has the 
expected results for constant J conditions. The results in 
Table-2 quantify the decrease in jet penetration with the 
increase in the number of jets as shown in Figure- 13. Note 
that the streamlines depicted in Figure-13 indicate almost lull 
turning of the jet by x/R=l and these results correlate well 
with the temperature and equivalence ratio distribution plots 
indicating minor developments in the mixing flow field 
beyond x/R=l. The temperature and equivalence ratio 
profiles however indicate that enhanced mixing is obtained 
with the increase in the number of orifices for a given value 
of J as reported previously by Oechsle, Mongia, and Holde- 
man. ( 1992). The results obtained with the area weighted, 
mass flow weighted, and histogram parameters B(+) and B(-) 
indicate this same conclusion. The difference in the histo- 
grams between 8 and 12 orifices/row at similar J. MR and DR 
conditions are shown in Figure- 14. Note that the histograms 
differ significantly in regions near f equl i (Figure- 14c and 
14d). The reacting flow configurations (Fig- 14b and 14d) 
show slightly less difference in mixing between 8 and 12 
holes/row (with cumulative Avol values of about 50%) as 
compared to the non-reacting configurations (Fig- 14a. and 
14c with Avol values in excess of 60%). 

As the number of holes increases, the jet wake volume 
decreases and therefore prevents the hot mainstream gas from 
entering and occupying this volume and reacting with cold 
jet stream which is recirculated in this wake. Both reacting 
and non-reacting results appear to be sensitive to the size of 
this jet trailing wake region as mentioned previously in this 
paper. The best mixing configuration from the 8 cases 
reported in this paper is the configuration number 7 shown in 
Figures-7g, 8g, and 1 lc. Note that the optimum mixer for a 
nominal J=25 is 9 holes/row as calculated in Eq-9 (Holde- 
man. 1993). where n is the number of orifices and C=2.5. 
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Effect of the rich-zone equivalence ratio on mixing 

The effect of the increase of the rich-zone equivalence ratio 
(tJtRZ) from 1.35 to 1.80 is shown in Figures-5b, 5d. 6b. and 
6d (temperature profiles) and Figures-7b. 7d. 8b. and 8d 
(equivalence ratio profiles). Significant differences in the 
temperature flow field are noted with the increase in 
rich-zone equivalence ratio. With a =1.35. the flow 
entering the mixing section need only mix slightly with the 
jet air to allow stoichiometric f/a to increase the temperature 
due to the expected chemical kinetics. However, the <J) RZ 
= 1.80 configurations require a longer time to allow mixing of 
the rich mixture entering the mixing section, and the resi- 
dence time of the mixing section limits the amount of further 
reaction occurring in the mixer. Therefore the temperature 
profiles are different as compared to the ()> RZ =1.35 configura- 
tions. This effect can be somewhat beneficial in the design 
of an RQL combustor since the <f» RZ appears to limit the 
location of the reaction between the rich mainstream flow 
and dilution jet mixing, thus allowing the designer to delay 
the reaction process after adequate mixing has been obtained 
in the mixing section. The difference in the amount of 
reaction in the mixing section is clearly shown in Figure-5f 
and 5h for (Jirz of 1.35 and 1.80 consecutively. In addition to 
the delayed reaction produced by increased in <j> RZ . the 
mixing flowfield is slightly enhanced (see the results of 
configurations 1 and 3. 2 and 4. 5 and 7. and 6 and 8). Note 
that the volume histograms cannot be compared for in this 
case since the f equi j is a function of the <)>R 2 • 

Conclusions 

1) The parameter (0 based on the non-reacting temperature 
profiles does not compare with the reacting flow normalized 
temperature solutions due to the additional reaction occurring 
in the mixing zone. The (tj parameter based on temperature 
and equivalence ratio however shows similar flow field 
solutions for non-reacting and reacting Hows (respectively) 
since both are conserved scalar quantities. 

2) The jet interaction between opposing jets becomes more 
significant for non-reacting Hows calculations as compared 
to reacting (low results however, jet penetration itself appears 
not to have been significantly affected. The core of the jet 
appears to diffuse less and has less overall mixing in the 
reacting flow solutions as compared to the non-reacting 
cases. 

3) All three statistical parameters, area weighted, mass flow 
weighted, and histogram half widths correlate well for the 
selection of the best mixing configuration (configuration #7). 
Significant differences are shown in the histogram shapes 
between reacting and non-reacting flows that are not shown 
using simple planar statistical parameters. This indicates that 
additional study is necessary to evaluate the use of planar 
averages in characterizing the mixing flow phenomena, and 
in particular, extrapolating the expected NO x formation in a 


mixing section of an RQL combustor based on results valued 
at x/R=l. 

4) Increase in the number of orifices decreases the overall 
volume of the wake behind each jet and therefore signifi- 
cantly affects the reaction in this area. The wake region 
behind the jet is an optimum location for flame stability due 
to low velocity recirculation areas where mainstream rich 
fuel mixture is mixed with relatively cold jet air to enhance 
reaction. 

5) An increase in rich-zone equivalence ratio (4 >rz) produces 
an effective reaction delay in the mixing section: therefore, 
the selection of the appropriate tf)^ allows the combustor 
designer to delay this mixing zone reaction sufficiently to 
allow adequate mixing to occur. 
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CASE 

Units 

i 

2 1 3 

4 

5 

6 

7 8 

OVERALL 










.1 (jet/main) 


25.0 

25.0 

25.0 

25.0 

25.0 

25.0 

25.0 

25.0 

MR (jet.main) 


1.83 

1.83 

2.67 

2.67 

1.83 

1.83 

2.67 

2.67 

DR (jet/main) 


2.578 

2.584 

2.300 

2.231 

2.578 

2.584 

2.300 

2.231 

VR (jet/main) 


3.114 

3.110 

3-297 

3.347 

3.114 

3.110 

3.297 

3-347 

Ar (jet/main) 


0.228 

1 0.228 

0.352 

0.358 

0.228 

0.228 

0-352 

0-358 

Phi (rich zone) 


0.00 

1J5 

0.00 

1.80 

0.00 

1-35 

0.00 

1.80 

Phi (lean zone) 


0.00 

0.45 

0.00 

0.45 

0.00 

0.45 

0.00 

0.45 

DP/P 

% 

2.91 

3.08 

3.00 

3.10 

3.00 

3.00 

3.00 

3.10 











MAINSTREAM 










P main 

atm. 

1.0 

13.4 

1.0 

13.4 

1.0 

13.4 

1.0 

13.4 

T main 

K 

773 

2-527 

690 

2,182 

773 

2,527 

690 

2.182 

V main 

m/sec 

22.7 

42-3 

2U 

39.4 

23.1 

41.7 

21.8 

39.4 

M main 

kg/sec 

0.052 

0.399 

0.056 

0.431 

0.045 

0-393 

0.056 

0.430 

Mixer Diameter 

m 

0.080 

0.080 

0.080 

0.080 

0.080 

0.080 

0.080 

0.080 











JET 


Circular 

Circular 

Circular 

Circular 

Circular 

Circular 

Circular 

Circular 

P jet 

atm. 

1.0 

13.4 

1.0 

13.4 

1.0 

13.4 

1.0 

13.4 

T jet 

K 

300 

978 

300 

978 

300 

978 

300 

978 

V jet 

m/sec 

70.8 

131.7 

71.9 

131.9 

71.9 

129.8 

71.9 

131.9 

M jet 

kg/sec 

0.096 

0.731 

0.150 

1.150 

0.083 

0.720 

0.150 

1.149 

ACd jet/row 

m 2 

1.15e-03 

1.14e-03 

1.77e-03 

lJI0e-03 

1.15e-03 

1.14e-03 

1.77e-03 

1.80e-03 

Orifice Cd 


0.80 

0.81 

0310 

0.82 

0.80 

0.80 

0.80 

0.82 

Number of orifices 


8 

8 

8 

8 

12 

12 

12 

12 

Blockage 


0.430 

0.430 

0.534 

0.538 

0.527 

0.526 

0.654 

0.659 

Orifice T.E. x/R 


0.338 

(U38 

0.420 

0.423 

0.276 

0.276 

0-343 

0J45 

Flow Condition 


NR 

REACT 

NR 

REACT 

NR 

REACT 

NR 

REACT 


Table-1. Overall description of the operating conditions and geometric configurations 


Conf 

Jet Penetration (a) x/R =0.625 

Jet Penetration <K x/R = 1.0 

Jet Penetration (ri x/R= 2.0 

w/R 

Penetration of jet core / R 

w/R 

Penetration of jet core / R 

w/R 

Penetration of jet core / R 

1 

0.425 

0.7 

035 

0.75 

0325 

0.75 

2 

0.4 

0.675 

035 

0.725 

0325 

0.725 

3 

0.475 

0.625 

0.45 

0.675 

0375 

0.7 

4 

0.475 

0.575 

0.425 

0.625 

035 

0.675 

5 

0325 

0.575 

0325 

0.625 

03 

0.65 

6 

035 

0.575 

0325 

0.625 

03 

0.65 

7 

0.425 

0.55 

0.4 

0.6 

0375 

0.6 

8 

0.425 

0.5 

0.425 

035 

0375 

0.575 


Table-2. Jet penetration development through the center of the jet in the circumferential direction 


Con fie 

condition 

# holes 

PHI RZ i MR 

A MIX 

A HOT 

ACOI.D 

MMIX 

MHOT 

MCOLD 

B(-) 

B(+) 

I 

NR 

8 

0 1.8.1 

0.21 

0.187 

0.096 

0.223 

0.195 

0.108 

0.171 

0.266 

2 

REACT 

8 

135 1.83 

0321 

0.191 

0.111 

0.223 

0.195 

0.108 

0.182 

0333 

3 

NR 

8 

0 2.67 

0.173 

0.154 

0.079 

0.184 

0.162 

0.085 

0.148 

0.255 

4 

REACT 

8 

1.8 | 2.67 

0.201 

0.175 

0.1 

0-209 

0.178 

0.109 

0.179 

0328 

5 

NR 

12 

0 

1.83 

0.164 

0.144 

0.079 

0.165 

0.138 

0.091 

0.157 

0.207 

6 

REACT 

12 

135 

1.83 

0.171 

0.144 

0.093 

0.171 

0.132 

0.108 

0.166 

0.244 

7 

NR 

12 

«_J 

2.67 

0.144 

0.13 

0.062 

0.139 

0.118 

0.074 

0.136 

0.169 

8 

REACT 

12 

1.8 

2.67 

0.166 

0.142 

0.082 

0.164 

0.131 

0.098 

0.159 

0.288 


Table-3. Summary of the mixing non-uniformity results for all the analyzed configurations 


NOTE: REACT - Reacting flow , NR = Non-reacting flow. PHI RZ = rich-zone equivalence ratio 
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Figure 1. Depiction of (a.) geometric configuration of the mixing section 
and (b.) typical computational grid. 
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RICH ZONE EQUIVALENCE RATIO AND 
MASS FLOW RATIO VARIATION 


Figure-2. Matrix of analyzed configurations 



Figure-3. Jet parameter definition 
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Figure-4. Graphical definition of B(+) and B(-) 
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1.5 

1.0 

'o. 


Configuration # I 
MR=I.83 y Non-reacting 


o.o 

- 0.5 
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x/R 
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Configuration # 3 
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Configuration # 7 
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MR- 2.67, Reacting 
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Figure-5. Normalized Temperature distribution through the center of the jet 

(a., b., c.y amt d. ) 8 round orifices / row jet = 0 

(e., f, g., amt //, ) 12 round orifices / row mainstream = t (non-react flows) 

stoichiometry = / (react flows) 
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a. Configuration # 1 
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Non-reacting 



0.0 0.5 1.0 
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c. Configuration # 3 
MR=2.67 
Non-reacting 



b. Configuration # 2 
MR=1.83 
Reacting 
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e. Configuration # 5 
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f. Configuration # 6 
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Reacting 



g. Configuration # 7 

MR-2.67 
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h. Configuration # 8 
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Figure-6. Normalized temperature distribution in the radial-tangential plane at x/R- 1 

(a.. I)., c., and d. ) S round orifices / row jet = 0 

(e.,f, g., and li. ) 12 round orifices / row mainstream = / (non-reacting flows) 

stoichiometry = / (reacting flows) 
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Figure-7 . Normalized Temperature distribution (/) {non-reacting flow) compared to the normalized 
equivalence ratio distribution (f) (reacting flow)— (plane through the center of the jet) 

(a., b„ c., and d. ) H round orifices / row f =0 at T n or 

(e., f, g., and It. ) 12 round orifices / row f- / at T tmin or 
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a. Configuration # 1 
MR= 1.83 
Non-reacting 
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c. Configuration # 3 
MR- 2. 67 
Non-reacting 



b. Configuration # 2 
MR= 1.83 


Reacting 
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li. Configuration # 8 
MR= 2.67 
Reacting 
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Figure-8. Normalized temperature (j) for the non-reacting flow compared to the normalized 
equivalence ratio (f) for the reacting flow for radial-tangential plane at x/R= I 


(a., b., c., and d. ) 8 round orifices / row 
( e.,f, g., and h. ) 12 round orifices / row 


f=0 at T jtl or 

f ^ 1 atT main 0t 
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Mixing non-uniformity ^ Mixing non-uniformity 


Area weighted Parameters 
results at x/R=l 


Mass flow weighted Parameters 
results at x/R=l 



Configuration number 


2 3 4 5 6 

Configuration number 


► 



1 2 3 4 5 6 

Configuration number 


s 


Figure-9. Mixing non-uniformity parameters (Table-3). 
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A Fractions for Conf (1-2) ft. A Fractions for Conf (3-4) 
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fl. A (7 - 2), M/?= 7.53, 5 holes/row b. A(3 - 4), MR=2.67, 8 holes/row 

c. A(5 - 6), MR= 1.83, 12 holes/row d. A(7 - 8), MR=2.67, 12 holes/row 

Figure-12. Difference between the reacting and non-react configurations for the same geometry and operating conditions 
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Figure-13. Jet penetration as seen through the center of the jet 
(a., b. y c and d. ) 8 round orifices / row 
(e.,f, g., and h. ) 12 rowid orifices / row 
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A Fractions for Conf (1-5) b. A Fractions for Conf (2 - 6) 
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a. A(I - 5), MR= 1.83, non-reacting b. A (2 - 6), \1R= 1.83, reacting 
c. A (3 - 7), MR=2.67, non-reacting d. A (4 - 8), MR=2.67, reacting 

Figare-14. Effect of the variation of the number of holes/row on mixing. 
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